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Abstract 

We study the role of massive perturbers (MPs) in deflecting stars and binaries to almost radial ("loss-cone") or- 
bits, where they pass near the central massive black hole (MBH), interact with it at periapse, and are ultimately 
destroyed. MPs dominate dynamical relaxation when the ratio of the 2nd moments of the MP and star mass 
distributions, /12 = N p {M^ j N± (M 2 ), satisfies (j,2 3> 1. We compile the MP mass function from published 
observations, and show that MPs in the nucleus of the Galaxy (mainly giant molecular clouds), and plausibly in 
late type galaxies generally, have 10 2 < \i% < 10 s . MPs thus shorten the relaxation timescale by 10 2 ~ 7 relative 
to 2-body relaxation by stars alone. We show this increases by 10 1-3 the rate of Zarge-periapse interactions with 
the MBH, where loss-cone refilling by stellar 2-body relaxation is inefficient. We extend the Fokker-Planck 
loss-cone formalism to approximately account for relaxation by rare encounters with MPs. We show that bi- 
nary stars-MBH exchanges driven by MPs can explain the origin of the young main sequence B stars that are 
observed very near the Galactic MBH, and increase by orders of magnitude the ejection rate of hyper- velocity 
stars. In contrast, the rate of sma/Z-periapse interactions of single stars with the MBH, such as tidal disruption, 
is only increased by a factor of a few. We suggest that MP-driven relaxation plays an important role in the 
3-body exchange capture of single stars on very tight orbits around the MBH. These captured stars may later be 
disrupted by the MBH via tidal orbital decay or direct scattering into the loss cone; captured compact objects 
may inspiral into the MBH by the emission of gravitational waves from zero-eccentricity orbits. 
Subject headings: black hole physics — clusters — galaxies: nuclei — stars: kinematics — giant molecular 
clouds 



1. INTRODUCTION 

There is compelling evidence that massive black holes 
(MBHs) he in the cen ters o f all galaxies (iFerrarese & Merrittl 
1200(1 iGebhardt et alJl2003UShields et alJl2003l). including in 
the ce nter of our Galaxy ( Eisenhau er et al.ll2005t iGhez et all 
120051) . The MBH affects the dy namics and evolution o f the 
galaxy's center as a whole (e.g. Bahca ll & Wolfl ll976) and 
it also strongly affects individual stars or binaries that ap- 
proach it. Such close encounters, which may be extremely 
energetic, or involve non-gravitational interactions, or post- 
Newtonia n effects, have be en the focus of many studies (see 
review by Alexander 2005). These processes include the de- 
struction of stars by the MBH, either by falling whole through 
the event hor izon, or b y being first tidally disrupted and then 
accreted (e.g.lReesll 988): tidal scattering of stars on the MBH 
J Alexander & Liviol20 01 ): the capture and gradual inspiral of 
stars into the MBH, accompanied by t he emission of gravita- 
tional waves or by ti dal heating (e.g. Alexander & Hopman 
120031 1 Alexander & Morri sl2003l) : or dynamical exchange in- 
teractions in which incoming stars or binaries energetically 
eject a star tightly bound to the MBH and are captured in 
its place very near the MBH (e.g. lAlexander & Livioll2004l 
Gould & Ouill enl2003l) . 

The interest in such processes is driven by their possible 
implications for the growth of MBHs, for the orbital de- 
cay of a MBH binary, for the detection of MBHs, for grav- 
itational wave (GW) astronomy, as well as by observations 
of unusual stellar phenomena in our Galaxy, e.g. the puz- 
zling young population of B-star very near the Galactic MBH 
jEisenhauer et aljf2005l) . or the hyper-velocity B stars at the 
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edge of the Galaxy dBrown et alJl2005l iFuentes et al1l2006t 
Brown et al. 2006a), possibly ejected by 3-body interactions 
of a binaries with the MBH (iHillsl 19911) . 

Here we focus on close encounters with the MBH whose 
ultimate outcome ("event") is the elimination of the incom- 
ing object from the system, whether on the short infall (dy- 
namical) time, if the event is prompt (e.g. tidal disruption or 
3-body exchange between a binary and the MBH), or on the 
longer inspiral time, if the event progresses via orbital decay 
(e.g. GW emission or tidal capture and heating). Such pro- 
cesses are effective only when the incoming object follows an 
almost zero angular momentum ("loss-cone") orbit with peri- 
apse closer to the MBH than some small distance q. To reach 
the MBH, or to decay to a short period orbit, both the infall 
and inspiral time s must be much shorter than the system's re- 
laxation time t r ( Alexander & Hopman 2003). The fraction 
of stars initially on loss-cone orbits is very small and they are 
rapidly eliminated. Subsequently, the close encounter event 
rate is set by the dynamical processes that refill the loss-cone. 

The loss-cone formalism u sed for estimating th e event 
rate (iFrank & Reesl Il976t iLightman & Shapiro! 119771 
ICohn&Kulsrudlll978l) usually assumes that the system is 
isolated and that the refilling process is 2-body relaxation. 
This typically leads to a low event rate, set by the long 2-body 
relaxation time. 

Two-body relaxation, which is inherent to stellar systems, 
ensures a minimal loss-cone refilling rate. Other, more effi- 
cient but less general refilling mechanisms were also studied 
with the aim of explaining various open quest ions (e.g. the 
stalling problem of MBH binary coalescence. IBerczik et alJ 
(2005); M erritt & Wand (120051): IBerczik et alJ (120061), or 
MBH feeding, Zhao et alJ2002UMiralda-Escude & Kol lmeier 
2005) or in the hope that they may lead to significantly higher 
event rates for close encounter processes. These mechanisms 
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includ e chaotic orbits in tr i axial p o tentials tNorma n & Silkl 
1 19831 Iderhard & Binnevl Il985t iMerritt & Pood 120041 
Hollev-Bockelmann & Sigurdsson 2006) (the presence of a 
MBH may howev er destr oy the triaxiality near the center; 
Merritt & Ouinlan 1998; iHollev- Bockelm ann et all 120021: 
Sellwood 2002); increased fraction of low angular momentum 
orbits in non-spherical potentials (Magorrian & Tremaine 
1 19991 iBerczik et all 120061) : accelerated resonant relaxation 
of angular momentum near the M BH whe re the orbits are 
Keplerian ilRauch & Tremainelll996|: iRauch & rngailslll998t 
Hopm an"& AlexandeJ l2006at iLevid 120061) : perturbations 
by a massive accretion d isk or an intermediate mass black 
hole (IMBH) companio n (iPolnarev & Reesll 19941 IZhao et al.1 
120021 iLevin et al.|l2005l) . Most of these mechanisms require 
special circumstances to work (e.g. specific asymmetries 
in the potential), or are short-lived (e.g. the IMBH will 
eventually coalesce with the MBH). 

Here we explore another possibility, which is more likely 
to apply generally: accelerated relaxation and enhanced 
rates of close encounters driven by massive perturbers 
(MPs) . Efficient relaxat i on by MPs w as first suggested by 
Spitzer & Schwarzschild (1951, 1953) to explain stellar ve- 
locities in the galactic disk. MPs remain an important com- 
ponent in modern models of galac t ic disk heating (see e.g . 
| VillumsedfT9^fT^lLacevlll98ltlJenkins & Binnevlll990t 
Hanninen & Flvnn 2002 and references therein). A similar 
mechanism was suggested to ex plain the spatial diffu sion of 
stars in the inner Galactic bulge (Kim & Morris 2001). In ad- 
dition to dynamical heating, efficient relaxation by MPs was 
suggested as a mechanism for loss cone replenishment and 
relaxation, both in the context of scatt ering of Oort cloud 
comets to the Sun dHillsl 1198 ll lBailevlfT983h an d the scat- 
tering o f stars to a MBH in a galactic nucleus ( Zh ao et all 
2002). IZhao et al. (2002) sugges ted MPs as a mechanism 
for est ablishing the M./cr relation (Ferrares e & Merriti2000t 
Gebh ardt et al.l l2000) by fast accretion of stars and dark mat- 
ter. They also noted the possibility of increased tidal disrup- 
tion flares and accelerated MBH binary coalescence due to 
MPs. In this study we investigate in detail the dynamical 
implications of loss-cone refilling by MPs. We evaluate its 
effects on the different modes of close interactions with the 
MBH, i n particular 3-body exchanges, which were not consid- 
ered by Zh ao et all (12002). and apply our results to the Galac- 
tic Center (GC), where observations indicate that dynamical 
relaxation is very likely dominated by MPs. 

This paper is organized as follows. In ^2] we present the 
main concepts and procedures of our calculations. The ob- 
servational data and theoretical predictions about MPs in the 
inner ~ 100 pc of the GC are reviewed in |3] In section |4] 
we explore the implications of relaxation by MPs for various 
types of interactions with the MBH. We summarize our results 

2. LOSS-CONE REFILLING 

In addition to stars, galaxies contain persistent dense struc- 
tures 3 such as molecular clouds, open clusters and globular 
clusters with masses up to 10 4 -10 7 M Q . Such structures can 
perturb stellar orbits around the MBH much faster than 2- 
body stellar relaxation (hereafter "stellar relaxation"), pro- 
vided they are numerous enough. This condition can be 
quantified by considering a test star randomly scattered by 

3 Structures that persist at least as long as the local galactic dynamical time 
and are substantially denser than the ambient stellar mass distribution. 



perturbers with masses in the interval (M p , M p + dM p ) and 
number density (dNp/dM p )dMp, approaching it with rela- 
tive velocity v on orbits with impact parameters in the inter- 
val (b, b+db). The minimal impact parameter still consistent 
with a small angle deflection is 6 m j n = GM p /v 2 (the capture 
radius), where v is of the order of the local velocity dispersion 
a. Defining B = b/b min > 1, the encounter rate is then 

( d 2 r \ JBrJl /< 



\dMpdb 



dM p db 



dN p 
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The total rate is obtained by integrating over all MP masses 
and over all impact parameters between &,„;„ and 6 max . Here 
we are interested in perturbations in the specific angular mo- 
mentum J of a star relative to the central MBH, and so 
b max ~ r, the radial distance of the star from the center. MPs 
with substantially larger impact parameters are much less ef- 
ficient because their effect on the MBH-star pair is tidal rather 
than direct. 

The relaxation rate due to all MPs at all impact parameters 
is then 



dM^db 



= dbJdM p ( 

logA$/dM p (^W, (2) 



where log A = log(6 max /6 m i n ) is the Coulomb logarithm 
(here the dependence of log A and v on M p is assumed to 
be negligible). For stars, typically logA> 10; the omission of 
large angle scattering by encounters with b < 6 m ; n is thus jus- 
tified because it introduces only a relatively small logarithmic 
correction. This formulation of the relaxatio n tim e is equiv- 
alent to its conventional definition (Spitzer 1987) as the time 
for a change of order unity in v 2 by diffusion in phase space 
due to scattering, t r ~ v 2 /D(v 2 ), where D(v 2 ) is the diffu- 
sion coefficient. 
If the stars and MPs have distinct mass scales with typical 

1 /2 

number densities AT* and N p and rms masses (M 2 y ,z and 

(M 2 ) 1/2 {(M 2 ) = j M 2 (dN/dM)dM/N), then MPs dom- 
inate if the ratio of the 2nd moments of the MP and star mass 
distributions, ^ 2 = N p (M 2 ) / AT* (M 2 ), satisfies [i 2 > 1 
(note that for a continuous mass spectrum, this condition is 
equivalent to -d log N/d log M < 2). 

As discussed in detail in |3j the central ~ 100 pc of the 
Galactic Center (GC) contain 10 8 — 10 9 solar masses in stars, 
and about 10 6 — 10 s solar masses in MPs such as GMCs 
or open clusters of masses 10 3 - 10 8 Mq (|Oka et alJl2001t 
[ Figer et all2002U2004UVollmer et all2003tlGusten"&P hilipp 
120041 IBorissov a et al.l2005l) . An order of magnitude estimate 
indicates that MPs in the GC can reduce the relaxation time 
by several orders of magnitude, 
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Note that /i 2 does not include possible modifications in the 
value of log A for MPs due to their much larger size, which 
may decrease this ratio by O(10). This estimate is borne 
by more detailed calculations (Fig. [2 and table |2}, us- 
ing the formal definition t r — v 2 /D(v 2 ^) with M*p* — > 
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J (dN p /dM p )M 2 dM p (e.g. iBinnev & Tremaindll987l Eqs. 
8-69 to 8-70). A similar result is indicated by simulat ions of 
spatia l diffusion of stars in the central 100 pc l lKim & Morris 
2001). 

2.1. Non-coherent loss-cone refilling 

The Fokker-Pl anck approach to the loss-cone problem (e.g. 
ICohn&Kulsrudl [l978) assumes that the effects of multiple 
small perturbations on the orbit of a test star dominate over the 
rarer strong close encounters (6 max /6 m i n 3> 1), and that the 
cumulative effect can be described as diffusion in phase space. 
The change in the angular momentum of the test star then 
grows non-coherently, A Joe y/t. The change over one orbital 
period P is Jd — J c (E)y / P/t r , where J c = \/2(ip — E)r is 
the maximal (circular) angular momentum for a stellar orbit of 
specific relative energy E~—v 2 /2 + ip(r), and if>= —<f) is the 
negative of the gravitational potential, so that E > for bound 
orbits. The magnitude of Jr> relative to the J-magnitude of 
the loss-cone, 

J lc ~y/2GM.q, (4) 

determines the mode of loss-cone refilling. The relative vol- 
ume of phase space occupied by the loss-cone, Jf c /J^{E), 
increases with E (decreases with r) while P decreases. Near 
the MBH (high E) J D -C J ic , stars diffuse slowly into the 
loss-cone, and are promptly destroyed over an orbital pe- 
riod, leaving the loss-cone always nearly empty. In this 
empty loss-cone regime, the loss-cone is relatively large, but 
the refilling rate is set by the long relaxation timescale (e.g. 
Lightman & Shapiro 1977), 
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where N+ (E) is the stellar number density per energy interval. 

Far from the MBH (low E) Jd ^> J; c , stars diffuse across 
the loss-cone many times over one orbit, and the loss cone 
is always nearly full. In this full loss-cone regime the refill- 
ing rate is set by the short orbital time, but the loss cone is 
relatively small, 



4 N*jE) 
J*(E) P(E) 



(6) 



Note that here and elsewhere we make the simplifying ap- 
proximation that the period is a function of energy only, which 
is true only for motion in a Keplerian potential. 

The total contribution to loss-cone refilling is dominated 
by stars with energies near the critical energy E c (equiva- 
lenfly, critical typical radius r c ) separating the two regimes 
jLiehtman & ShapirolH977l see 5g). Within r c (E > E c ), 
an object, once deflected into the loss cone, can avoid being 
scattered out of it before reaching the MBH 4 . The empty and 
full loss-cone regimes of infall processes can be interpolated 

4 In the case of inspiral, E c is determined by the condition Jjj = J; c over 
the inspiral time, rather than the much shorter orbital period, which results in 
a much smaller r c than for direct infall. Inspiraling stars with E > E c can 
avoid being scattered directly into the MBH before completing the orbital 
decay. There is no contribution to inspiral events from regions outside r c 
(E < E c ), since the probability of an object to remain on its low- J trajectory 
over the many orbital periods required to complete the inspiral, is vanishingly 
small. 



to give a general approximate expression for the differen- 
tial ev ent rate for these non-coherent encounters (e.g. lYoung l 
fl977h . 



dT 
dE 



with 



3 2 {E) 



j 2 (E) N*(E) 
J?(E) P(E) ' 

J D (E) 
log(J c (E)/J lc ) 



•Ji, 



(7) 



(8) 



where j is the loss-cone limited angular momentum change 
per orbit, which expresses the fact that the loss-cone can at 
most be completely filled during one orbit. 

2.2. Coherent loss-cone refilling 

The loss-cone formalism can be generalized to deal with 
MPs in an approximate manner with only few modifications. 
The capture radius for MPs may be smaller than their size. 
Since the MP mass profile is centrally concentrated, we adopt 
the modified definition 



b min = mm(0.1R p ,GM p /v 2 ) 



(9) 



where v 2 = GM(< r)/r. This results in log A ~ 6 for MPs 
, less than the typical value for relaxation by stars. Never- 
theless, the error introduced by neglecting encounters with 
b< 6 m in is still not very large because penetrating encounters 
are much less efficient. However, the assumption of multiple 
non-coherent encounters with MPs over one orbital period is 
not necessarily justified because of their small number den- 
sity. 

To address this, we modify the treatment of the empty loss- 
cone regime (the contribution to the event rate from regions 
where the loss-cone is already filled by stellar relaxation can 
not be increased by MPs, see ^4}. We define rare encounters 
as those with impact parameters b<b\, where bi is defined by 



bi 



P / db{dT/db) = l 



(10) 



The differential rate is estimated simply by (dT/db) — 
N p v2nb. When P f*™* d6(dT/d6) > 1, with b max = r, all 
encounters with b > b\ are defined as frequent encounters 
that occur more than once per orbit, and add non-coherently 5 . 

Note that even when P J h bmax dfc(dr/d6) > 1 for all MPs, per- 
turbations by rare, very massive MPs may still occur less than 
once per orbit. Our treatment is approximate. A complete sta- 
tistical treatment of this situation lies beyond the scope of this 
study. 

When the typical number of encounters per orbit is less than 
one, the fractional contributions from different individual en- 
counters, 8 J, should be averaged coherently (AJoci), subject 
to the limit that each encounter can at most fill the loss-cone. 
The loss-cone limited change in angular momentum per orbit 
due to rare encounters is therefore 



3%(E) = 



bl dT 
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In contrast, frequent uncorrected collisions add up non- 
coherently (A J oc v£)> an d it is only their final value that 



In the marginal cases of P J b max d6(dr/d6) 



< 1 or 



P / bi max db(dr/db) < 1, all encounters 



are considered rare. 
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is limited by the loss-cone (individual steps 6 J may exceed 
Jic, but can then be partially cancelled by opposite steps dur- 
ing the same orbit). The loss-cone limited change in angular 
momentum per orbit due to frequent encounters is therefore 



j F (E) =min 



1 



log(J c /J /c ) 



P 



*db^5J 2 ,J? c 



(12) 



The total loss-cone limited angular momentum change per or- 
bit is then approximated by 

i 2 = min(j|+i|,jf c ) , (13) 

and the differential event rate is calculated by Eq. |7}, 
AT/AE = [f(E)/J*(E)] N*(E)/P(E). 

The contribution of rare encounters is evaluated in the 
impulse approximation by setting 5J ~ GM p r/bv in Eq. 
fu\ . We find that the this contribution by GC MPs (50 
is generally small. Frequent encounters are the regime usu- 
ally assumed in the Fokker-Planck treatment of the loss- 
cone problem (e.g. Lightman & Shap irolll977l) . To evalu- 
ate the contribution of frequent encounters, we do not cal- 
culate 5 J directly, but instead calculate the sub-expression 

I = P J^ x d6(dr/d6)(5J 2 in Eq. <0 in terms of the fr- 
ay eraged diffusion coefficient D(v 2 ), after averaging over the 
orbit between the periapse r p and apoapse r a and averaging 
over the perturber mass function (this i s ess entially equivalent 
to the definition of Jd in terms of t r , 32.H . 
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The assumptions involved i n the last approximate term 
(Magorrian & Tremaine 1999) are that the star is on a nearly 
radial orbit (v r — ► v, r p — > 0, r a — ► 2r) and that D(v 2 ) (the 
diffusion coefficient of the transverse velocity relative to the 
MBH) can be approximated by D(Av? L ) (the diffusion coef- 
ficient of the transvers e velocity relative to the stel lar velocity 
v), given explicitly by (Binnev & Tremaine ( 1987), Eq. 8-68) 

„ 87rG 2 (diVp/dM p )M B 2 lnA / v \ 
D(A V 1) = \ * K (-^ , (15) 

where K{x) = erf(ir)(l — 1 /2a; 2 ) + exp(— x 2 ) /y/wx and 
where a spatially homogeneous distribution of MPs with a 
Maxwellian velocity distribution of rms ID velocity a was 
assumed. 

To summarize, the event rates are calculated as follows. For 
each perturber model (table |2j, we integrate over the stellar 
distribution (N± = 1.2 x 10 6 (r/0.4pc) -2 , for r > 0.4pc and 
M* = 1 Mq) in terms of r, using 2V* to deri ve the appropri- 
ate d ensity of perturbed objects (single stars 34.1l or binaries 
^4.21 . At each r we calculate 6 min (Eq. |9}, bi (Eq. [lOj, j R 
(Eq. II lb and jp (Eq. I12> . The integral / (Eq. 1141 is evaluated 
by taking v 2 — > GM(< r)/r and correcting approximately 
for the difference relative to the exact calculation. We use j 
(Eq. 1 1 3i to calculate the differential event rate AT/Ar (Eq. 
Q, with the ansatz E -> GM(< r) /2a where a = (4/5)r 
is an effective semi-major axis, which is motivated by the 
fact that for a Keplerian isothermal eccentricity distribution, 
(r) = a(l + (e 2 ) /2) = (5/4)a. The total event rate is cal- 
culated by carrying the integration over r in the region where 
MPs exist, between tmp and r out (§j3j. 
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FIG. 1 . — Relaxation time as function of distance from the MBH, for stars 
(solid line) and for each of the 4 MP models separately, as listed in table |2| 
clusters (dashed-dotted lines), GMCs (dashed lines). The discontinuities are 
artifacts of the assumed sharp spatial cutoffs on the MP distributions. 2-body 
stellar processes dominate close to the MBH, where no MPs are observed to 
exist. However, at larger distances massive clumps (at 1.5 < r < 5 pes) and 
GMCs (at 5 < r < 100) are much more important. 



3. MASSIVE PERTURBERS IN THE GALACTIC CENTER 

MPs can dominate relaxation only when they are massive 
enough to compensate for their small space densities. Here 
we consider only MPs with masses M p > 10 2 Mq. Such 
MPs could be molecular clouds of different masses, in par- 
ticular giant molecular clouds (GMCs), open or globular stel- 
lar clusters, and perhaps also IMBHs. As discussed below, 
observations of the Galaxy reveal enough MPs to dominate 
relaxation in the central 100 pc. We adopt here a conserva- 
tive approach, and include in our modeling only those MPs 
that are directly observed in the Galaxy, namely GMCs and 
young clusters. We discuss briefly theoretical predictions for 
two other classes of MPs: dynamically evolved "submerged" 
clusters and IMBHs, that could well be common in galactic 
centers and contribute to efficient relaxation. 

The dynamically dominant MPs are GMCs. Emission line 
surveys of the central ~ 100 pc reveal ~ 100 GMCs with 
estimated masses in the range 10 4 -5 x 10 7 M Pl and sizes 
of Rr, ~ few pc jMivaza ki & Tsuboi 2000; Oka et alJl200lt 
Giisten & Philipp 2004). We selected individual, reliably 
identified GMCs in the central 0.7° of the Galaxy (^ 100 p c 
of the GC), from the sample observed by Oka et al. (2001). 
Figure |2 shows the empirical GMC mass function using their 
virial mass estimates as an upper mass limit and adopting a 
lower limit 10 times smaller, following Mivazaki & Tsuboi 
(2000) who found that LTE mass estimates are typically an 
order of magnitude lower than the virial ones. Note that 
the more r ecent GMC CO 1-0 molecular line observations by 
lOka et akl (1200 ll) . which we use here, indicate a more mas- 
sive and flatter mass function than that derived for their ear- 
lier C S1-0 molecular line observations (Mivazaki & Tsuboi 
2000). This is probably due to the higher sensitivity of the 
CO 1-0 line to lower-density molecular gas (M. Tsuboi, priv. 
comm.). 

Inside the inner 5 pc of the GC (a volume smaller or compa- 
rable to that of a GMC or a stellar cluster), the most massive 
local structures are the molecular gas clumps observed in the 
circumnuclear gaseous disk (CND) and its associated spiral- 
like structures (Genzeletal. 1985; Christopher et al. 2005). 
Their size are ~ 0.25 pc and their masses are estimated to 
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be in the range 10 3 -10 5 Mq, where the lower estimates are 
based on the assumption of optically thin HCN(l-O) line emis- 
sion and the upper estimates are based on the optically thick 
assumption, which also coincides with the virial estimates. 

It is possible to obtain a model-independent estimate of 
the effect of the observed MPs on the relaxation time by 
the directly derived value (j,?, bs = J2i M p,i/ N * M *' which 
is listed in table [2] The observed GMC masses show that 
/i° bs - 2 x 10 6 ~ 8 on the 100 pc scale, and ^2 bs ~ 6 x lO 1 " 3 on 
the 5 pc scale, a clear indication that MPs dominate the relax- 
ation on all relevant lengthscales. where MPs exist. For the 
purpose of our numeric calculations below, it is convenient 
to describe the differential mass function analytically. Here 
we adopt a power-law dN p /dM p oc M^ 13 parameterization 
(or a lognormal probability distribution function (LNPDF) in 
the appropriate cases dN p /dM p oc LNPDF(fi, a), where/i 
is the log mean and cris the log standard deviation). Figure 
12 shows that the pow er law d istibution is a goo d fit for the 
GMCs mass function (Mivazaki & Tsuboi 2000), and for the 
lower estimates of the clusters and gas clumps. For the up- 
per estimates of the clusters and the molecular gas clumps the 
MFs were found to be better fit by a lognormal distribution, 
which we used in these cases. It should be emphasized that 
our results and conclusions are determined primarily by the 
large values of (/,%, and not by the detailed form of the mass 
function. We repeated our calculations with several alterna- 
tive distributions and found qualitatively similar results. The 
assumed high mass cutoff of the MF is important, as it deter- 
mines the magnitude of Figure |2] shows that our models 
do not extrapolate beyond the maximum observed MP masses 
(and even fall below them for GMCs). 

We obtained best fit power indices (3 and LNPDF indices 
/land a for the lower and upper estimates. The cumulative 
mass functions and best fits are shown in Fig. and listed 
in tabled We find (3 = 1.2 for both the lower mass range 
(1.4 x 10 4 < M p < 5 x 10 6 Mq) and the upper GMC mass 
range (1.4 x 10 5 < M p < 5 X 10 7 M Q ). The clump mass 
function has (3 ~ 1.1 for the lower mass range (2.4 x 10 2 < 
M p < 1.1 x 10 4 M o ) and fx = 10.04, a = 0.65 (/3~ 1.7 for 
the best power law fit) for the upper mass range (3.6 x 10 3 < 
M p < 1.35 x 10 5 M Q ) (Fig. Hand table[TJ. The space density 
of such clumps falls rapidly inside the inner ~ 1.5. 

Stellar clusters are another class of MPs, which are of 
minor dynamical significance in this context. About 10 
young stellar clusters with masses in the range 10 2 -10 5 Mq 



fled (pisreretal.1119991 12002 


Maillard et al. 2004; Figer et al. 


2004; Borissovaet all 12005 


). Again, we 


St the lower and 



upper mass estimates of these clusters (with mass ranges 
3 x 10 2 - 1.3 x 10 4 M Q and 4.5 x 10 2 - 7 x 10 4 Mq, respec- 
tively) with a power-law mass function of a LNPDF, and find 
j3 ~ 1.3 and fj, ~ 8.65, a = 1.1 ((3 = 1.9 for the best power 
law fit; see Fig. [2] and table 0. It is interesting to note in 
passing that both the GMCs and gas clumps and the clusters 
have bery similar mass functions distributions, but the clus- 
ters have a lower mass range, as might be expected if these 
GMCs are the progenitors of young GC clusters, similar to the 
relatio n observed be tween galactic disk clusters and GMCs 
JLada & Ladal2003l) . 

Based on the current observations of the 9 confirmed clus- 
ters in the GC, we find that they are dynamically insignifi- 
cant compared to the GMCs. Note however that their con- 
tribution to the relaxation is comparable to that of the stars, 




10 10 

Mass (M a ) 

FIG. 2. — Cumulative mass functions of observed massive perturbers in the 
GC (symbols) with best fit power-law differential distributions, dN/ dM oc 
M~ ''(dashed lines) or lognormal distribution with log mean [land log stan- 
dard deviation cr (solid lines). Top panel: lower (circles) and upper (squares) 
estimates for the masses of observed molecular clumps i Christopher et al. 
2005 ) (left) and GMCs IQka et alJ200!h (right) in the GC . The best fit power- 
law indices for the GMCs are (3 up = fti ow = 1.2 and for the clumps f3 up = 
1.1 and Pi ow = 1.7. For the upper estimates of the clumps masses a better 
fit is found by a log normal distribution wit h fi = 8.2, a = 0.65. Bottom 
panel: likewise for GC clusters i Figer et al. 1999; Figer 2004; Borissova et al. 
I200S) with /3 up = 1.3 and 0i ow = 1.9, and fi = 10, a = 1.1 for the upper 
masses estimate. 



even when the lower mass estimate is assumed. When the 
upper mass estimate is assumed, clusters shorten the relax- 
ation time by a factor of 10 (Fig. [2 and table |5Jl. Clusters 
could play a larger role in relaxation if there are many more 
yet undi scovered ones in the GC. B ased on dynamical sim- 
ulations, Porteg ies Zwart et"al] (1200 6) suggest that there may 
be O(100) evolved young clusters in the region, undetected 
against t he field stars beca use of their low surface density. 
However. lFiger et all (120021) argue that these massive clusters, 
it they exist, should have been observed. Two young G C clus- 
ter can didates were recently discovered by Bo rissova et alJ 
(2006), one of which may be more massive than 10 4 M©. 
It is thus possible that the GC harbors some additional un- 
detected massiv e clusters, although p robably not as many 
as suggested by Portegies Zwart et al. (2006). Such young 
cluste rs may grow IMBHs by runaway s tellar collisions (e.g. 
Portegies Zwart ~alJ2004tlFreitag et all2004|) . which would 
then sink to the center by dynamical friction, dra gging with 
them the cluster core. Portegies Zwart et al. (2006) estimate 
that many IMBHs may have migrated to the central 10 pc of 
GC in this way, however there is as yet no observational evi- 
dence supporting this idea. 

The MBH's dynamical influence extends up to a radius rn 
containing 0(M,) mass in stars, which in the GC falls in the 
range 2-4 pc, depending on the operative definition of rn and 
the uncertainties in the value of tf, and of the de nsity of the 
surrounding stellar cluster (see I Alexander! 120051 sec 3.1.2). 
Thus most Galactic MPs, and in particular the massive ones, 
lie outside r^. Table Q summarizes the estimated properties 
of the observed MPs . 

The observed MP species vary in their spatial distributions 
and mass functions, which are not smooth or regular, and their 
distribution constantly changes as their orbits decay by dy- 
namical friction and new MPs are formed and destroyed. It is 
thus likely that the presently observed MPs are but one real- 
ization of a much smoother and regular (well-mixed) under- 
lying distribution. Consequently, we construct several simpli- 
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Table 1 

Abundances of observed massive perturbers in the Galactic Center 

MP type r"(pc) N p M V (M & ) (3 (n, a) (Afg) 1/2 (M ) B. v (pc) References 

Observed GMCs < 100 -100 10 4 - 10 s L2 3 X 10° -3 X 10' 5 Oka et al. (2001); Glisten & Philipp (2004) 

Observed clusters < 100 -10 10 2 - 10 5 1.3 (8.2, 1.1) 4.8 X 10 3 - 2.4 X 10 4 1 Figer et al. (1999. 2002) : Figer (2004) 

M.ull.lld Ol l.:iH)4 l. lidllVMA,! CI ill, ijmi^ i 

Observed clumps 1.5 - 3 -25 10 2 - 10 5 1.1 - 1.7 (10, 0.65) 3.7 X 10 3 - 4.1 xlO 4 0.25 Genzel et al. (1985); Christopher et al. (2005) 
"Projected distance range enclosing observed MPs. 



Table 2 
Massive perturbers models 



Model 


r (pc) a 


N p » 


M p (M Q ) 






GMC1 


5-100 


100 


10 b -5 x 10 ' 


1.2 






1.5-5 


30 


3 x 10 3 - 10 5 


1.1 


10.04, 0.65 


GMC2 


5-100 


100 


10 4 -5 x 10 6 


1.2 






1.5-5 


30 


7 x 10 2 -10 4 


1.7 




Clusters 1 


1.5-100 


10 


5 x 10 2 -7 x 10 4 


1.3 




Clusters2 


1.5-100 


10 


3 x 10 2 -10 4 


1.9 


8.16, 1.1 


Stars 


5-100 


2xl0 8 


1 






Stars 


1.5-5 


8xl0 6 


1 







"Distance range enclosing observed MPs. 
b N p (r) oc r -2 in the given range is assumed. 

c M2 b3 = T.i M p j/ N * M ^ where M p,i is the observed MP's mass. 



fied MP models for our numeric calculations (table|2j that are 
based on the observed properties of the MPs (|4), and our best 
fits for their mass functions. The simplifications involve the 
following assumptions, (i) A smooth, spherically symmetric 
MP number density distribution, N p (r) oc r~ 2 between the 
cutoffs tmp < r < r out (r out does not play an important role, 
see below), with a random velocity field, (ii) A single or bro- 
ken power-law MP mass functions, dN p /dM p oc M~@ (or a 
lognormal distribution in the appropriate cases), (iii) A single 
mass stellar population of 1 M@ stars with a number density 
distribution N*(r) oc r~ 2 outside the inner 1.5 pc (i.e. a con- 
stant MP to star ratio), (iv) Mutually exclusive perturber types 
(i.e., a single type of perturber is assumed to dominate relax- 
ation, as indicated by the detailed calculations presented in 
Fig. 

The five MP models are detailed in table|2] Stars, Clusters 1, 
Clusters2, GMC1 and GMC2, represent respectively the case 
of relaxation by stars only, by heavy and light stellar clusters, 
and by heavy GMCs and light GMCs. Table|2]lists /i 2 , the 
ratio of the 2nd moment of the various MP mass distributions 
to that of the stars. 

4. MASSIVE PERTURBER-DRIVEN INTERACTIONS WITH A MBH 

The maximal differential loss-cone refilling rate, which is 
also the close encounters event rate, dF/dE, is reached when 
relaxation is efficient enough to completely refill the loss cone 
during one orbit (Eq. 0. Further decrease in the relaxation 
time does not affect the event rate at that energy. MPs can 
therefore increase the differential event rate over that pre- 
dicted by stellar relaxation, only at high enough energies, 
E > E c (equivalently, small enough typical radii, r < r c ), 
where slow stellar relaxation fails to refill the empty loss- 
cone. The extent of the empty loss-cone region increases with 
the maximal periapse q, which in turn depends on the close 
encounter process of interest. For example, the tidal disrup- 
tion of an object of mass M and size R occurs when q < r t , 



Binary 
Stars 


MPs 

E Relaxation 




No MPs 


J v 




E 


E ', Stellar % 
■ Relaxation % 


F 


Single Stars 


_EJ F 




E 


I >v F 





n " mr Radius (log scale) c 

FIG. 3. — A schematic representation of the local contribution to the tidal 
disruption rate for a single component MP model (e.g. GMC2 without the 
central clumps, so that rjjp = 5 pc). Top: the disruption rate of binaries 
(a constant binary fraction is assumed here for simplicity) due to stellar re- 
laxation (bottom dashed line) and MPs (top dashed line). Bottom: the dis- 
ruption rate of single stars for both perturber models. Empty and full loss- 
cone regimes are denoted by "E" and "F", respectively. The initial orbital 
energy of the disrupted objects is expressed by the corresponding radius of 
a circular orbit, r. For r < r^, stellar density and velocity distributions of 
N* (r) oc r~ 3 / 2 and cr(r) oc r -1 / 2 are assumed; for r>Tf l , N* (r) oc r~ 2 
and cr(r) = const. The transition from the empty to full loss-cone regime 
(for stellar relaxation) occurs at the critical radius < r^p for single stars 
and at ^> tmp for binaries. 



the tidal disruption radius, 

r t ~R(M./M) 1/3 



(16) 



This approximate disruption criterion applies both for sin- 
gle stars (M — M*, R = i?*) and for binaries, where M is 
the combined mass of the binary components and R is the 
binary's semi-major axis, a. Stellar radii are usually much 
smaller than typical binary separations, but stellar masses are 
only ~ 2 times smaller than binary masses. Binaries are there- 
fore disrupted on larger scales than single stars. In the GC this 
translates to an empty (stellar relaxation) loss-cone region ex- 
tending out to r s c ~ 3 pc for single stars and out to r c > 100 
pc for binaries. In the GC tmp < r s c -C r h c , and so MPs 
are expected to increase the binary disruption rate by orders 
of magnitude, but increase the single star disruption rate only 
by a small factor. This is depicted qualitatively in figure Q, 
which shows the local rates (dr / d log r) for disruption of sin- 
gle stars and binaries due to stellar relaxation or relaxation by 
a simplified one component MP model. 

Figure also shows that the MP-induced disruption rate is 
dominated by binaries originating near the inner cutoff tmp 
(in the following discussion the initial orbital energy of the 
disrupted objects is expressed by the corresponding radius of 
a circular orbit, r ~ GM(< r)/2E). This is qualitatively dif- 
ferent from the usual case of stellar tidal disruption induced by 
stellar relaxation, which mainly occurs inside r s c and is dom- 
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inated by stars ori ginating near the outer limi t min(r;?, r^); 
usually r s c ~ rn (Li ghtman & Shapiro! 1 1977t> . The differ- 
ence can be understood by considering the r-dependence of 
dr /d log r . Neglecting logarithmic terms, the empty and full 
local loss-cone rates are, respectively (Eqs.|5j|6j 



dr e _ 7V*(<r) 
dlogr t r {r) 

dlogr 



M, 



M. + M±(<r) 



q\ N+{<r) 
P(r) 



(I) 



(17) 



(18) 



where N*(< r) is the number of stars enclosed within r. 

Inside Th, the potential is dominated by the MBH (M*(< 
r)/M. -C 1) , A^oc r~ ai wither < 2 for most dy namical sce- 
narios (e.g. iBahcall & Woll 19771 lYoun^ fl980) and the ve- 
locity dispersion is Keplerian, a oc r -1 / 2 . The orbital period 
scales as Pocr 3 / 2 and t r ~ (M./M*) 2 P/[(log/V*)iV*]. The 
local disruption rates are then dT e /dlogr 



r<r h 



OC r 



9/2-2qi 



oc r 1 / 2 ai . For plausible values of 
, increases with r whereas Tf de- 



lr>r h 
„-2 



cxr - 2 / 2 ~ 
For «2 ~ 2 



and dr//dlogr| r<n> 
1/2 < ai < 9/4, T, 
creases with r, so the rate is dominated by stars near r u 
l lLightman & Shapiroll 19771) . Outside r/j, the stellar distribu- 
tion is usually assumed to be near isothermal, oc r~" 2 with 
ai ~ 2 and a velocity dispersion er ~ const, and the potential 
is dominated by the stars (M, <C M*(< r) cx r 3 ~" 2 ). The 
orbital period scales as Pcxr" 2 / 2 and t r ^N±P/ logiV*. The 
local disruption rates are then dr e /dlogr 
r- 1 and dT/ /d log r | r>r& cx r -i-<W2 _ 
both rates increase toward small radii. Since for most Galactic 
MP types rMP > Th, the disruption rate is dominated by stars 
near rMP- For example, when the loss-cone is empty, ~ 50% 
of the total rate is due to MPs at r < 2rMp; when the loss-cone 
is full, ~ 75% of the total rate is due to MPs at r < 2r M p. 

4.1. Interactions with single stars 

GMCs, gas clumps and clusters in the GC are abundant only 
beyond the central rMP ~ 1.5 pc, whereas the empty loss- 
cone regime for tidal disruption of single stars extends only 
out to r s c ~ 3 pc. For inspiral pro cesses such as GW emission, 
r r is ~ 100 times smaller still (Hopm an & Alexanderll2005l 
2006b). The effect of such MPs on close encounter events 
involving single stars is thus suppressed (weaker tidal effects 
by MPs at r > r s c a re not considered here). This is contrary 
to the suggestion of Zhao et al. (2002), who assumed that the 
effect of MPs fully extends to the empty loss-cone regime. 
We find that the enhancement of MPs over stellar relaxation 
to the single stars disruption rate is small, less than a factor 
of 3, and is due to stars scattered by gas clumps in the small 
empty-loss cone region between rup ~ 1.5 pc and r s c ~ 3 pc. 
A possible exception to this conclusion is the hy pothesized 
population of IMBHs (Portegies Zwart et al. 2006), not mod- 
eled here, whose distribution co uld extend to the inner pc (e.g. 
the IMBH candidat e in IRS 13.|Maillard etld~ll2004l but see 
iSchodel et al .120051 andlPaumard et all 

4.2. Interactions with stellar binaries 

The empty loss cone region for binary-MBH interactions 
extends out to > 100 pc because of their large tidal radius. 
On these large scales MPs are abundant enough to dominate 
the relaxatio n processes . Here we focus on 3-body ex change 
interactions dHillslll 99 U 119921 lYu & Tremainel2 003 ), which 
lead to the disruption of the binary, the energetic ejection of 



one star, and the capture of the other star on a close orbit 
around the MBH. 

Various phenomena associated with such excha n ge in- 
teractions were suggeste d and e xplored. IHillsl dl988|) 
and later lYu & Tremainel J2Q03I). iGualandris et al"l (2005), 
Ginsburg & Loebl (120061) and iBromlev et alJ ll200 6). studied 
hyper-velocity stars ejected from the GC fol lowing tidal dis- 
ruption by the MBH. iGould & Ouillenl ( 120031) suggested this 
mechanism to explai n the origin o f the young stars near the 
Galactic MBH. Mill er et all (12005) proposed that compact ob- 
jects captured following a binary disruption event will eventu- 
ally be sources of GWs from zero-eccentricity orbits, in con- 
trast to high-eccentricity sour ces typical of single star inspiral 
(Hopman & Alexander 2005). 

The event rates estimated by these authors vary substan- 
tially. IHillsl (11988) assumed a full loss-cone and a fraction 
/bin = 0.02 of the stars in binaries with small enough semi- 
major axis to produce a high-velocity star (a < 0.1 AU), and 
derived a 3-body exc hange rate of ~ 10~ 3 (/bin/0.02) yr _1 . 
lYu & Tremainel (120031) took into account the empty loss cone 
regime, and argued for a higher fraction of relevant bina- 
ries (/bin = 0.04 for binaries with a < 0.3 AU that can 
survive 0.8 pc from the MBH), thereby obtaining a rate of 
~ 2.5 x 10 _6 (/bi n /0.04) yr _1 , 3 orders of magnitude smaller 
than that estimated by Hills. These calculations assumed the 
same binary separation for all binaries and a constant binary 
fraction at all distances from the MBH ( two possibilities were 
considered, a = 0.3 AU and a = 0.03 AU). 

The binary fraction and typical binary semi-major axis de- 
pend on the binary mass, and on the rate at which binaries 
evaporate by encounters with other stars. This depends in turn 
on the stellar densities and velocities, and therefore on the dis- 
tance from the MBH. Here we take these factors into account 
and estimate in detail the 3-body exchange rate for MP-driven 
relaxation. The rate is proportional to the binary fraction in 
the population, which is the product of the poorly-known bi- 
nary IMF in the GC and the survival probability against binary 
evaporation. 

We assume for simplicity equal mass binaries, A/bin = 
2M*, since the observations indicate that moder ate mass ra- 
tios dominate the binary pop ulation ( iDuauennov & Mayor 
1199 It iKobulnickv et alJl2006). The evaporation timescale at 
distance r from the center for a binary of semi-major axis a 
compo sed of two equal mass stars of mass M* and lifetime t* 
is (e.g. lBinnev & Tremainel 1 9871) 



t e vap(M*,a,r) = 



bin 



o-(r) 



(19) 



(M*) 160Fp(r)GalnA bin 

The Coulomb factor for binary evaporation, Abi n = 
aa 2 1 4G (MJ) , expresses the fact that the binary is only af- 
fected by close perturbations at distances smaller than ~ 
a/2. The MPs considered here are extended objects (table|2j 
and therefore do not affect the binary evaporation timescale 
(IMBH MPs could be a possible exception). Although bi- 
nary evaporation is a stochastic process and the actual time 
to evaporation differs from binary to binary, we expect a 
small scatter in the evaporation rate, At~^ ap / t~^ p <C 1, be- 
cause evaporation is a gradual process caused by numerous 
weak encounters. Evaporation is thus better approximated as 
a fixed limit on the binary lifetime, rather than as a Pois- 
son process (where At~^ ap / t^} ap = 1). The maximal bi- 
nary lifetime is then r = min([tjj, t t (M ± ), t cvap (M+, a, r)], 
where tn is the Hubble time, taken here to be the age of the 
Galaxy. It is well established that the central 100-200 pc of 



8 



Perets, Hopman and Alexander 



the GC have undergon e continuous star formation over the 
lifetim e of the Galaxy ( Sera bvn & MorrisH l996: Fig er et all 
l200l . Assuming a constant star formation rate over time tn, 
the differential binary distribution at time t is dN^/ da\ t = 
/bin d//da| 1\ mm(t, t), where d//da| is the normalized 
initial semi-major axis distribution, which can be observed in 
binaries in low-density environments where £ cva p — > oo, and 
r* is the single star formation rate, which is normalized to the 
observed present day stellar density by setting t = tn and tak- 
ing i evap — >oo for singles, so that T+ = N ic / min(£jj, The 
present-day binary semi-major axis distribution is therefore 



dN hii 



da 



=Ain(M*; 



d/ 
da 



N*(r) x 



min < 1, 



£cva P (A^,a,r) 
min[ijr,t*(M*)] 



(20) 



The capture probability and the semi-major axis distribu- 
tion of captured stars were es timated by simulations (Hills 
1 1 99 ll fl992t I Yu & Tremainel20 03). Numeric experiments in- 
dicate that between 0.5-1.0 of the binarie s tha t approach the 
MBH within the tidal radius r t (a) (Eq. 1 161 are disrupted. 
Here we adopt a disruption efficiency of 0.75. The harmonic 
mean semi-major axis for 3-body ex changes with equal mass 
binaries was found to be (lHillsll991l) 



0.56 



M m 



2/3 



0.56 



/ Mm 



\M hh 



1/3 



n, (21) 



where a is the semi-major axis of the infalling binary and 
a\ that of the captured star (the MBH-star "binary"). Most 
values of a\ fall within a factor 2 of the mean. This rela- 
tion maps the semi-major axis distribution of the infalling 
binaries to that of the captured stars: the harder the bina- 
ries, the more tightly bound the captured stars. The veloc- 
ity at infinity of the ejected star (neglecting the Galactic po- 
tential) is u| H = 2 1 / 2 (GM hh Ja)(M./M hin ) 1 / 3 oc M^Ja 
(an eq ual mass binary with periapse at r t is assumed; Hills 
1988). The harder the binary, the higher is wbh- The pe- 
riapse of the captured star is at r *, and therefore its eccen- 
tricity is very high dHillsl I199U [1992; Mi ller et all 120051) . 
e = 1 - r t /ai ~ 1 - 1.8(M bin /M.) 1 / 3 > 0.98 for values 
typical of the GC. 

We now consider separately the implications of 3-body ex- 
change interactions of the MBH with old (t* > tn) binaries 
and massive young (t ± < 5 x 10 7 yr) binaries. 

4.2.1. Low mass binaries 

The properties of binaries in the inner GC are at present 
poorly determined. The period distribution of Solar neigh- 
borhood binaries can be approximated by a log normal dis- 
tributi on with a median per iod of 180 years (a ~ 40 AU) 
(DuQu emiov & Mavorlll99il) . The tot al binary fra ction of 
these binaries is estimated at /bin ^0.3 (Lada 2006). Adopt- 
ing these values for the GC, the total binary disruption rate 
by the MBH can then be calculated, as described in ^2] by 
integrating diVbin/da (Eqs. I20i over the binary a distribution 
and over the power-law stellar de nsity distribution of the GC 
up to 100 pc JGenzel et alJl2003l) . Table (0 lists the capture 
rates for the different perturber models, assuming a typical old 
equal-mass binary of Mbj n = 2 M@. 

The old, low-mass binary disruption rate we derive for stel- 
lar relaxation alone is ~ 5 x 10~'yr -1 , ~ 5 times lower, 
but still in broad agreement with the result of lYu & Tremaind 



(2003). Their rate is somewhat higher because they assumed 
a constant binary fraction and a constant semi-major axis for 
all binaries, even close to the MBH, where these assumptions 
no longer hold. 

MPs increase the binary disruption and high-velocity star 
ejection rates by factors of ~ 10 1-3 and effectively acceler- 
ate stellar migration to the center. This can modify the stel- 
lar distribution close to the MBH by introducing a "source 
term" to the stellar current into the MBH. Low-mass stars 
are at present too faint to be directly observed in the GC. 
However, such a source term may have observable conse- 
quences since it can increase the event rate of single star pro- 
cesses such as tidal disruption, tidal heating and GW emis- 
sion from compact objects, in particular from compact objects 
on zero-eccentricity orbits ( Mille r~t al J 120051) (in contrast, 
GW from inspiralin g single stars o ccur on high-eccentricity 
orbits, Hopman & Alexander 2005). We calculated numer- 
ical solutions of the Fokker-Planck equation for the stellar 
distribution around the MBH with a captured stars source 
term. These preliminary investigations (Perets, Hopman & 
Alexander, in prep.) confirm that the total accumulated mass 
of captured stars does not exceed the dyna mical cons traints 
on th e extended mass around the MBH ( Mouawad et al. 
120051) . bec ause 2-body relaxation and likely also resonant 
relaxation (Rau ch & Tremainel 119961 iHopman & Alexanderl 
2006a) scatters enough of them into the MBH or to wider or- 
bits. 



4.2.2. Young massive binaries 

MPs may be implicated in the puzzling presence of a 
cluster of main sequence B-stars (4 < M* < 15 Mq) in 
the inner ~ l" (~ 0.04 pc) of the GC. This so-called "S- 
cluster" is spatially, kinematically and spectroscopically dis- 
tinct from the young, more massive stars observed farther 
out, on the ~ 0.05-0.5 pc scale, which are thought to have 
formed fro m the gravitational fragmen t ation of one or two 
gas disks (Levin & Beloborodov 2003; Genzel et a fl l2003l 
iMilosavlievic & Loebll2004l iPaumard et al.ll2006l) . There is 
however still no satisfactory explanation for the existence of 
the seemingly normal, young massive main sequence stars of 
the S- cluster, so close to a MBH (see review of proposed mo d- 
els bv lAlexanderr2 005: also a recent model by Levin 2006). 

Here we revisit an idea proposed by [Gould & Ouillen 
J200l . that the S-stars were captured near the MBH by 
3-body exchange interactions with infalling massive bina- 
ries. Originally, this exchange sc enario lacked a pl ausible 
source for the massive binaries. iGould & Ouiilenl specu- 
lated that they originated in an unusually dense and mas- 
sive young cluster on an almost radial infall trajectory, but 
concluded that such a finely-tuned scenario seems unlikely. 
Furthermore, a massive cluster is expected to leave a tidally 
strippe d tail of massive stars be yond the central 0.5 pc 
jKim et alJl2004t IGiirkan & Rasiol 120051) . which are not ob- 
served ( Paumard et al. 2006). Alternatively, it must con- 
tain an unusually massive central IMBH to hold it together 
again st the tidal field of the GC JHansen & Milosa vlievid 
2003). However, such a massive IMBH is well beyond what is 
predict ed by simulations o f IMBH formation by runaway col- 
lisions dGiirkan et alJl2004l IGiirkan & Rasiol2005l). or antici- 
pated by ex trapolating the Ml a relation ( iFerrarese & Merrittl 
l2TMlGebhardt et all200(]l) to clusters. 
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Table 3 

Total binary disruption rate and number of captured young stars 



Model 


Disruption rate (yr 1 ) 
r<0.04pc r<0.4pc 


Young Stars" 
r<0.04pc 


Young Stars" 
0.04 < r<0.4pc 


Young HVSs" 


GMC1 


1.1 x lCT* 3.2 x 1CT 4 


39.4 


5.2 


342 


GMC2 


2.8 x icr 5 1.8 x itr 4 


6 


1.4 


49 


Clusters 1 


6.4 x 1CT 7 9.7 x itr 7 


0.19 


0.004 


1.9 


Clusters2 


2.6 x icr 8 4 x itr 8 


0.01 


lO" 4 


0.1 


Stars 


3.4 x 1CT 7 5.3 x 10~ 7 


0.15 


0.003 


1.3 


Observed 


? ? 


10 - 35 c 


? 


43±31 d 



"Main sequence B stars with lifespan t < 5 X 10 Y yr . 

b Main sequence B stars with lifespan t < 4 X 10 8 yr. Notice that only 20 percents of these ejected stars could be observed in regions covered by current surveys. . 
c ~10 stars with derived a< 0.04 pc. >30 stars are observed in the area. 



Estimated from the observed 5 HVSs, at distances between 20-120 pc from the GC I Brown et al. 2006: 



MP-driven 3-body exchanges circumvent the problems of 
the cluster infall scenario by directly bringing massive field 
binaries to r t , without requiring massive clusters of unusual, 
perhaps even impossible properties. The ongoing star forma- 
tion in the central ~ 100 pc implies the presence of a large 
reservoir of massive stars there, which are indeed observed 
in the central few x 10 pc both in dense clusters and in the 
field JFieer et aljfl999t lKgerll2003t iMuno et al.l l2006). It is 
plausible that a high fraction of them are in binaries. 

We model the binary population of the GC in the S-stars 
mass range , 4 < M± < 15 M Q , by assuming equal mass bi- 
naries that follow the single star mass function with an ini- 
tial binary frac tion of fhin ~ 0.75, as observed elsewhere 
in the Galaxy (lLadall2006l iKobulnickv et al.ll2006l) . Because 
the stellar evolutionary lifespan of such stars is relatively 
short, massive binaries are essentially unaffected by dynami- 
cal evaporation. We assume star format ion at a constant rat e 
for 10 Gyr with a Miller-Scalo IMF ( Mille r & Scalolll979i) 
and u se a stellar population synthesis code ( Sternberg et al. 
2003) with the Geneva stellar evolution tracks ( Schall er et alJ 
1992) to estimate that the present day number fraction of stars 
in the S-star mass range is 3.5 x 10 -4 (and less than 0.01 
of that for M* > 15 M @ stars). Note that if star forma- 
tion in the GC is biased toward massive stars (Figerl l2003l 
IStolte et 10120051) . this estimate should be revised upward. 
We adopt the observed Solar neighborhood distribution of 
the semi-major axis of massive binaries, which is log-normal 
with (logo) = -0.7±0.6 AU (i.e. 63% of t he binaries with 
a = Q 2jjn i° AU; 91% with a = 0.2±^| AU) dGarmanv et all 
1 19801 iKobulnic kv et all20 06). Massive binaries are thus typ- 
ically harder than low-mass binaries, and will be tidally dis- 
rupted (Eq. I16> closer to the MBH and leave a more tightly 
bound captured star. 

We represent the massive binaries by one with equal mass 
stars in the mid-range of the S-stars masses, with A/bin = 
15 M and f*(7.5 M®) ~ 5 x 10 7 yr, and integrate over the 
stellar distribution and the binary a distribution as before, to 
obtain the rate of binary disruptions, T, the mean number of 
captured massive stars in steady state, AT* = Tt+, and their 
semi-major axis distribution (Ea.l21>. Table (0 compares the 
number of captured young stars in steady state, for the differ- 
ent MP models, on the r < 0.04 pc scale (the S-cluster) and 
0.04 < r < 0.4 pc scale (the stellar rings) with curre nt obser- 
vations (Eisenh auer et all2005UPaumard et al.l2006 j). 

The S-stars are found in th e central < 0.04 pc (iGhez et aTl 
120051 Eisenhau er et al.ll2005l) . If they were captured by bi- 
nary disruptions, they must have originated from massive bi- 



naries with a < 3.5 AU. This is consistent with semi-major 
axis distribution of massive binaries. The number of captured 
massive stars falls rapidly beyond 0.04 pc (table because 
wide massive binaries are rare. This capture model thus pro- 
vides a natural explanation for the central concentration of 
the S-cluster (Fig|4j. The absence of more massive stars in 
the S-cluster (M* > 15 A/ Q , spectral type O V) is a statistical 
reflection of their much smaller fraction in the binary popula- 
tion. Figure (0J and table (0 compare the cumulative semi- 
major axis distribution of captured B-stars, as predicted by the 
different MP models, with the total num ber of young stars ob- 
served in the inne r 0.04 pc (~ 35 stars. [ Eisenhauer et alJ200"5l 
IGhez et alJl2005l iPaumard et al.ll20o4 . Of these, only ~ 10 
have full orbital solutions (in particular a and e) at present. 
For the others we assume the ansatz that a is similar to the 
observed projected position. The numbers predicted by the 
GMC-dominated MP models are consistent with the obser- 
vations, unlike the stellar relaxation model that falls short by 
two orders of magnitude. 

The binary capture model predicts that captured stars have 
very high initial eccentricities. Most of the solved S-star 
orbits do have e > 0.9, but a couple have e ~ 0.3- 
0.4 JEisenhauer et alJ 120051) . Normal, non-coherent stel- 
lar relaxation is slow, even after taking into account the 
decrease in t r toward the ce nter due to mass segregation 
( Hopman & Alexander 2006b). It is unlikely that it could 
have decreased the eccentricity of these stars over their rel- 
atively short lifetime s. However, the much f aster process of 
resonant relaxation (iRauch & Tremainelll996l) may be effi- 
cient enough to randomize the eccentricity of a fraction of 
the stars, and could thus possibly expla in the much larger ob - 
served spread in eccentricities (Hopman & Alexander 2006a). 
Additional orbital solutions and a better estimate of the effi- 
ciency of resonant relaxation in the GC are required for more 
detailed comparisons between observations and the MP model 
predictions. 

4.2.3. Hyper-velocity stars 

Each captured star is associated with an ejected compan- 
ion, which in some cases is launched with a very high ve- 
locity. The one-to-one correspondence between the num- 
ber of captured S-stars and the number of early-type hyper- 
velocity stars (HVSs) is thus a generic prediction of binary 
capture models. The MP capture scenario specifically im- 
plies the continuous and isotropic ejection of both young 
and old HVSs from the GC. Recent observations of HVSs 
jBrown et all2005UHirsch et alJl2005l lEdelmann et al.ll2005t 
Brow n et al!2006albl) are consistent with a GC origin and fa- 
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FIG. 4. — Cumulative number of young B-stars in the GC as predicted 
by the MP models and by stellar two-body relaxation (listed in table 2). The 
ve rtical bar represents the total number of observed young stars inside 0.04 
pc lEisenhauer et al. 2005; S. Gillessen, priv. comm.). The dotted vertical 
line marks the approximate maximal distance in which captured B-stars are 
expected to exist. 



vor a steady state temporal distribution and an isotropic spatial 
distribution over a burst-like non-spherical distribution that is 
expected for HVSs triggered by the infall of a c luster (iLevinl 
l2005llHaardt et al.l2 006: Baumgar dtet alJ20 06). 

Two of the recently observed HVSs (v ~ 560- 
710 kms" 1 , Brown et al] 120051 I2006at Ententes et al.ll2006t 
Edelman net alJ 120051 iBrown et alJ 12006b ) were spectrally 
identified as late B V young massive stars (masses of 3-5 M©, 
MS lifespans of (l-4)xl0 8 yr and number fraction ~ 10 -3 in 
the population), implying a total popul ation of 43 ± 31 such 
hyper- velocity stars in the Galaxy (Br own et alJl2006ah . We 
model the parent binaries of these HVSs by equal mass bina- 
ries of M bin = 8 Mq and t*(4M Q ) = 2 x 10 8 yr The ejec- 
tion velocity was fou nd in numerical simulations dHillsll988t 
iBromlev et al.l2006T) to scale as 



«bh = 1776 kms 



0.1 AU 



-i/2 (Ml 



2Mr. 



1/3 



M 



BH 



3.7xl0 6 M(r 



1/6 

(22) 



To reproduce the high HVS velocities we consider binaries 
with a < 1 AU, which are tidally disrupted at r t < 3.7 x 
10~ 4 pc and eject a HVS wit h i>rh 900 k ms' 1 , the escape 
velocity from the bulge (Haardt et al. 2006). Taking only the 
GMCs into account, we predict that tens to hundreds such 
HVSs exist in the Galaxy, in agreement with the deduced HVS 
populations, whereas stellar relaxation predicts only 1.3 such 
stars (see table|3}. 

We note Brow neTaTl d200 6a) used the 10~ 5 yr -1 total rate 
of hyp er- velocity star ejection calculated by lYu & Tremaind 
( 2003) (including binaries of all stellar types, assuming stellar 
relaxation only and normalized to a fiducial 10% binary frac- 
tion) to estimate the number of late B V HVSs in a Salpeter 
IMF at 10-25. This theoretical prediction seems in rough 
agreement with the observations (and contradicts our much 
lower estimate of 1.3 HVSs). However, the rates of Yu & 
Tremaine are inapplicable here and lead to a significant over- 
estimate of the number of HVSs because their binary popula- 
tion model is not appropriate for massive binaries in the GC. 
On the one hand, the young binary population does not extend 
all the way to the center, as assumed by Yu & Tremaine for 
the general binary population (Following Paumar d'et alJ2006l 



who do not find any B V stars between 0.5-1 pc of the GC, 
we truncate the massive binary population inside 1.5 pc). The 
massive binary population in the GC is fewx 10 times smaller 
than implied by a simple scaling of the Yu & Tremaine gen- 
eral binary population. On the other hand, the binary fraction 
of young massive binaries is 70% rather than 10% 04.2.21 , 
We conclude that the agreement found by Brown et al. is 
accidental, and that binary disruption by stellar relaxation is 
insufficient to explain the number of observed HVSs, whereas 
MP-induced relaxation can reproduce the observations. 

5. SUMMARY 

Relaxation by MPs dominates relaxation by 2-body stel- 
lar interactions when the ratio between the 2nd mo- 
ments of their respective mass functions satisfies /12 = 
N p (Mf) I Ni, {Ml) > 1. We show that Galactic MPs (stel- 
lar clusters, GMCs and smaller molecular gas clumps that 
exist outside the inner few pc) dominate and accelerate re- 
laxation in the inner ~ 100 pc of the GC. This is plausibly 
the case in the centers of late-type galaxies in general. There 
is also evide nce for molecular gas in th e centers of late-type 
galaxies (e.g. Rupen 1997; Knapp 1999), which suggests that 
MPs may dominate relaxation there as well, and lead to the 
relaxation of the central regions of galactic bulges in general. 

Relaxation determines the rate at which stars and binaries 
are deflected to near radial (loss-cone) orbits that bring them 
closer to the MBH than some critical periapse q, where they 
undergo a strong destructive interaction with it. The size of q 
depends on the nature of the interaction of interest (e.g. tidal 
disruption, 3-body exchange). It is much larger for binaries 
than for single stars due to the binary's larger effective size. 

We extend loss-cone theory to approximately treat rare en- 
counters with MPs, and apply it to explore the implications of 
MPs on the rates of different types of close encounters. The 
rate reaches its maximum when loss-cone orbits are replen- 
ished by scattering within an orbital time (the full loss-cone 
regime). This is more easily achieved when the phase-space 
volume of the loss-cone is small, that is, when q is small. MPs 
thus affect only those processes with large q whose loss-cone 
is too large to be efficiently replenished by stellar encounters 
(the empty loss-cone regime). 

We show that MPs will not contribute much to the dis- 
ruption of single stars in the GC, whose loss-cone is effi- 
ciently replenished by stars outside the central ~ 2 pc (MPs 
may accelerate the consumption of stars by more massive 
MBHs, where q is significantly larger, or the capture of stars 
in accretion disks). However, MPs will enhance by fac- 
tors of 10-1000 the tidal disruption rate of infalling binaries, 
which result in the capture of one of the stars on a tight or- 
bit around the MBH, and the ejection at high velo city of the 
other star (iHillsll 1991lfl992t I Yu & Tremainell2003l) . The en- 
hancement of the event rates is dominated by the innermost 
MPs, and so the uncertainty in the MP distribution on the 
smallest scales dominates the uncertainties in the total event 
rate. Detailed observations of MPs in the inner GC allow 
us to robustly predict their effects in the Galaxy. We show 
that MP-induced disruptions of relatively rare massive bina- 
ries can naturally explain the puzzling presence of normal- 
app earing main sequence B stars in the central 0.04 pc of the 
GC ( Eisen hauer et alJ2005l) . and at the same time can account 
for the observed HVSs well on the ir way out of the Galaxy 
JBrown et all2005UE delmann et alJl2005UHirsch et al.ll2005t 
Brown et al. 2006a b; Bromlev et al. 2006). Tidal disruptions 
of the many more faint low-mass binaries can efficiently sup- 
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ply single stars on very eccentric tight orbits near the MBH. 
Such an increase in the numbers of stars in tight orbit near 
the MBH may increase the rates of single star processes such 
as tidal disruption and heating or GW emission from tightly 
bound compact objects (Miller et al.l2005l) . 

Finally, MP-induced interactions also have cosmological 
implications for the coalescence o f binary MBHs following 
galactic mergers JZhao et all2002t) . We suggest that MPs can 
accelerate the dynamical decay of binary MBHs by efficiently 
supplying stars for the slingshot mechanism, and thereby help 
solve the "last parsec" stalling problem. MP-driven loss-cone 
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refilling will operate even in the case of a spherical poten- 
tial, where other suggested mechanisms are inefficient, thus 
allowing MBHs to coalesce on the dynamical timescale of the 
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